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ABSTRACT 



Context. The halo of NGC 891 has been the subject of studies for more than a decade. One of its most striking features is the large 
asymmetry in Her emission. In this letter, we will take a quantitative look at this asymmetry at different wavelengths for the first time. 
Aims. We suggest that NGC 891 is intrinsically almost symmetric and the large asymmetry in Ha emission is mostly due to dust 
attenuation. We will quantify the additional optical depth needed to cause the observed asymmetry in this model. 
Methods. By comparing large strips on the North East side of the galaxy with strips covering the same area in the South West we can 
quantify and analyze the asymmetry in the different wavelengths. 

Results. From the 24 /im emission we find that the intrinsic asymmetry in star formation in NGC 891 is small i.e., ~ 30%. The 
additional asymmetry in Ha is modeled as additional symmetric dust attenuation which extends up to ~ 40" (1.9 kpc) above the plane 
of the galaxy with a mid-plane value of t=0.8 and a scale height of 0.5 kpc 
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1. Introduction 

NGC 891 is one of the best studied edge-on galax- 
ies in the nearby universe . Its close proximity (9.5 Mpc, 
; Ivan der Kruit & Searld [l98"Tl) and high inclination (89.8+0.5°, 
iKregel & van der Kruit] [2005) make NGC 891 an excellent can- 
didate for the study of the vertical structure of spiral galaxies. 

The z-structure of the emission from NGC 891 has been 
studied at many wa velengths over the past ~ 30 years. In bands 
such as the optical (van der Kruit & Searle, 1981; Xilouri s et all 
1998) and the HI (ISancisi & AllenL Il979h Is waters et all 119971) . 
NGC 891 is almost completely symmetric. The blue colors in 
the plane are brighter on the the North East side of the galaxy. 
This agrees with the idea of a trailing spiral structure with the 
dust behind the HII regions and the fact that the North East side 
is the approaching sid e as fi rst suggested by R. Sancisi (see, 
Ivan der Kruit & Searle ( 1981 ) a nd ref e rences therein). Th e first 
sensitive Ha studies (iDettmail Il990t iRand et all 1 19901) how- 
ever, observed a large asymmetry between the North-East (NE) 
and the South- West (SW) of NGC 891, extending up to large 
distances above the plane. Ha is a tracer of ionized gas, often 
linked to star formation (SF from here on) because generally 
only young, massive stars produce enough ionizing photons to 
ionize large volumes of neuttal gas. 

This asymmetry in He could be an effect of dust attenuation 
affecting this tracer in the same way as the blue colors in the 
plane mentioned earlier. Even high above the plane (|z|- 2 kpc) 
dust s tructures can be located and identified dHowk & Savagd 
2000) and it is natural to think that attenuation by this dust on 
the SW side of NGC 891 c an cause the a s ymme try. 

lAllen et all (1 19781) and iDahlem et all Jl994l) observed NGC 
891 with high resolution and high sensitivity in radio contin- 
uum emission, finding the asymmetry observed in Ha to be 
prese nt to a lesser extent in the distribution of radio contin- 
uum. [Dahlem et all (1 19941) suggest that this is the result of out- 



flow of relativistic electrons produced as a by-product of SF in 
the disk. This and the fact that the asymmetry in Ha extends 
well above the pla ne has led people to believe that it must be 
an intrinsic effect dDahlem et all 119941: iHowk & Savage! 120001: 
lKamphuisetaill2007h . " 

Recently a new S F tracer unaff e cted b y dust attenuation 
has become available. ICalzetti et all d2005l) show that 24 pm 
emission has a tight correlation with star formation. Although 
this correlation does not wor k in the same way for all galaxies 
dPerez-Gonzalez et all 120061) . the total Far Infrared Luminosity 
of NGC 891 (log(L(FIR) )= 10.27 L H ) is so cl ose to that of M51 
(log(L(FIR))= 10.42 L , ISanders et al.ll2003T) that this con-ela- 
tion can be used for NGC 891. Based on this new SF tracer we 
will provide for the first time a quantitative multi-wavelength 
analysis of the asymmetry in the halo of NGC 891. 



2. Data 

From the Spitzer Archive we retriev ed the public mosaics 
(PBCD) of NGC 891 at 3.6 pm (IRAC) (lAshbv et all lin preol 
iHolwerda et all lin prep .1) and at 24 pm (MIPS), with the archive 
tool Leopard. 

For the analysis of Ha primarily the image of IRand et ail 
(1990) is used. An integr ated map const r ucted from a Ha 
Fabry-Perot measurement (Kamphuis et all 120071) is used to 
rule out any systematic effects which might be present in the 
imaging. Since the cube is cal i brated on the image we will 
refer to the image of Ra nd et ail (l990) in the discussion unless 
mentioned otherwise. 

The Near and Far Ultra- Violet images of NGC 891 are 
retrieved from the GALEX public archive through the MAST 
websitdZl Unfortunately, the NUV is too crowded with fore- 
ground stars to give any useful information on diffuse emission 



http://archive.stsci.edu 
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Name 


(Central) Rest Wavelength 


Type of emission 


FWHM 


noise (<x) 




HI 


21 cm 


Line 


30" 


0.261 


MJy sr~' 


21 cm 


21 cm 


Continuum 


17" 


0.005 


MJy sr~' 


24 /urn 


24 /urn 


Continuum 


6.4" 


0.033 


MJy sr~' 


3.6 yum 


3.6 yum 


Continuum 


2.6" 


0.59 


MJy sr~' 


Hor (Rand et al.. 1990) 


6563 A 


Line 


1.2" 


7.8xl0- 18 


erg cirr 2 s _1 arcsec~ 2 


Ha (Kamohuis et al.. 2007) 


6563 A 


Line 


2" 


6.0xl0-' 8 


erg cm -2 s -1 arcsec -2 


NUV 


2275 A 


Continuum 


5.1" 


0.003 


MJy sr~' 


FUV 


1550 A 


Continuum 


4" 


0.002 


MJy sr~' 



Table 1. Summary of the characteristics of the different wavelength images. 
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Fig. 1. Upper Panel: Ha distribution of NGC 891 from iRand etaTI d 19901) in units of 10" 18 erg cm" 2 
to the 6.4"resolution of the 24 /mi. Overlaid in black contours is the 24 /mi emission at 3, 6, 12, 24 cr levels with cr=0.033 MJy 
sr _1 . Lower Panels: Intensity profiles parallel to the minor axis, averaged over 180"in the radial direction and 6.4"in the vertical 
direction. The profiles are normalized to the peak emission on the NE side. The panels show the NE (left panel) and SW (right 
panel) side of the galaxy. Solid blue line: Ho-, Dashed red line: 24 /mi. Note that the difference between NE and SW in 24 /mi is 
minimal. 



out of the plane of the galaxy. 

The HI map and 2 1 cm continuum image are from the most 
recent observations of N GC 891 with the Wes terbork Radio 
Synthesis Radio telescope dOosterloo et al.Ll2007l) . 

Except for blanking stars in all the optical and IR images, no 
further reduction was performed on the images. The resolution 
of the optical and IR images was determined by fitting gaussian 
point spread functions to 6 stars in each image. The main 
characteristics of the images are summed up in Table Q] 

Figure Q] (upper panel) shows the distribution of Ho- 
emission overlaid with contours of the 24 /mi emission. The 
lower panel of Fig. Q] shows the intensity profiles parallel to the 
minor axis, averaged over 180"in the radial direction and 6"in 
the vertical direction. The profiles are normalized to the peak 
emission on the NE side. The NE side of NGC 891 is shown 
in the left panel and SW side in the right panel. The solid line 
represent the Hff and the dashed line the 24 /mi. In Fig.Q]the 24 
//m profiles appear to be much narrower than the Hff profiles. 
However, at low levels these profiles have large extended wings. 
These wings are still well above the noise level of the image but 
are hardly seen in the vertical profiles due to the large dynamical 
range of the emission. 



3. Results 

We are interested in determining what fraction of the SF in 
NGC 891 is intrinsically asymmetric, and whether the remainder 
could be reasonably explained as extinction owing to dust at high 
|z| in the galaxy. We have averaged over large strips on either 
side of the galaxy to obtain the highest S/N possible. The strips 
have a width of 180"son the major axis on either side and start 
at a radial offset of 10"from the center so that effects from a 
central active source will not contaminate our results. We are 
interested in how the ratio between the emission in the NE and 
the SW changes with height in the different wavelengths. The 
resolution of the 24 /mi (i.e., ~ 6 ", see Table [TJ is the lowest 
resolution in our different bands with the exception of the 21 cm 
data. Therefore we have set the height of the strips equal to one 
resolution element in the 24 /mi image. The average intensity of 
such a strip on the NE is then divided by the average intensity of 
the strip on the SW with the same vertical offset. This way we 
obtain a vertical profile for the ratio between NE and SW with a 
resolution of 6". These vertical profiles are shown in the top half 
of Figure [2] 
Figure|2]shows us that: 

1. The emission in 24 /mi, 3.6 /mi and HI is almost symmetric. 

2. Ha and UV emission is asymmetric in the plane and declines 
to the same asymmetry as 24 //m, 3.6 /mi and HI above the 
plane. 
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Fig. 2. Upper Panels: Ratio betwe en the total emiss i on (see text) in the Nor th East and the South West. Left panel FUV (Galex), NUV 
(Galex); Left middle panel Ha dRand et all ll990t)Kamphuis et all 120071) : Right middle panel 24 /mi (Spitzer), 3.6 /mi (Spitzer); 
Right panel HI, 21 cm continuum dOosterloo et al.Ll2007l) . Al l panels show t wo do t-dashed lines for (Gray) a ratio of 1 indicating 
complete symmetry and (Black) a ratio of 1.3 as found by iDahlem et alj (1994) in the 6 cm radio continuum. Lower Panels: 
Additional t between t he spiral arms as derived from the Ha. Solid blue lines as from lRand et al.l dl990l) . Dashed red lines as from 
Kamphui s et al.l d2007l) . Left Panel: assuming complete symmetry. Right Panel: assuming an intrinsic asymmetry as observed in the 
24 /mi. 



We see that the ratio of asymmetry is very different between 
the wavelengths which are affected by dust attenuation (e.g., Ha 
and the UV) and the ones which are not (e.g., HI, 24 //m and 3.6 
/mi). Also, as we move away from the plane of the galaxy the 
asymmetry becomes less pronounced in the bands affected by 
dust attenuation. 

The 24 /mi and HI emission are not affected by dust 
attenuation and the ratio between the NE and SW is constant 
with height in these bands. The ratios in the plane of the 
galaxy are 1.24, 1.03 and 1.19 for the 24 /mi, HI and 21 
cm continuum respectively. The ratio found in the 24 /mi 
emission is very similar to t he ratio of 1.3 found in the radio 
continuum (A — 6 cm) by Dahlem et a D d!994h . Therefore, 
we can consider this ratio as the intrinsic asymmetry of SF 
in NGC 891. We model the rest as the effects of dust attenuation. 



4. Discussion 

Given the symmetry of the gas (HI), non-attenuated stellar 
(3.6 /mi) and SF (24 /mi) emission (Fig [2] right middle panel), 
all tracing different components of NGC 891, we assume that 
NGC 891 is intrinsically symmetric in the plane. Furthermore, 



we assume that the Ha above the plane is related to the SF in the 
plane. This means that the Ha above the plane is also symmetric 
and situated above the HII regions in the plane. We constructed 
a model based on the trailing spiral arm idea, which means that 
the HII regions are located in front of the dust lanes of the spiral 
arms. Then, because the SW is the receding side of the galaxy, 
on the NE we would be looking first at the HII regions and then 
the dust in the arm but on the S W side our line of sight to the HII 
regions passes through the arm first. Figure [3] shows the face-on 
view of our model of NGC 891. 

For the attenuating dust in the halo we assume an isotropic 
distribution (e.g., not specifically localized above the SF re- 
gions). It is already shown that the distribution of attenuating 
dust and ionized gas ar e unlikely to trace the same component 
dHowk & SavageL l2000). Furthermor e, on the SW side tendr ils, 
chimneys and shells are identified bv lHowk & Savage] d2000h . It 
is unlikely that these structures would be backlit by the stellar 
disks if they resided in the back of the galaxy. 

Therefore, in our model, the asymmetry in the Ha could 
easily be explained by dust in the halo. If the spiral arm on the 
SW side would be further away from us than the spiral arm on 
the NE, then the sight line through the dust toward the Ha on 
the South West of the galaxy would be much longer. From this 
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NE SW 



HII regions 




Fig. 3. Face-on view of our model for a symmetric NGC 891. 
The black S shape symbolizes the dust lane in the spiral arms 
and the blue areas the HII regions. The straight double headed 
arrow shows the difference in line of sight between the North 
East side and the South West side. The curved single headed 
arrow indicates the direction of rotation in NGC 891. 

idea we can easily calculate the additional attenuation between 
the two arms. The lower panel in Figure [2] shows the additional 
optical depth in the halo between the two spiral arms. r a dd is cal- 
culated by assuming that the Ha emission from the NE side of 
the galaxy is equal to the intrinsic emission on SW side and that 
all the differences are caused by dust attenuation. This is shown 
in the lower left p anel of Figu re [2] where the solid black line 
is from the image dRand et all |1990|) and the solid gray l ine is 
from the integrated moment map dKamphuis etal ., 2007]). The 
dust extends easily to 50"(2.3 kpc) above the plane. Even if we 
take the small asymmetry of the 24 //m in the plane into account 
(Fig 13 lower right panel) additional attenuation is present up to 
~40"(1.9kpc). 

Fitting a single exponential to the inner part (|z| < 20") of 
these optical depth profiles, we obtain a scale height for the dust 
of 0. 67 ±0. 10 kpc if the galaxy is perfectly symmetric and of 0.46 
+ 0.10 kpc if the intrinsic asymmetry is ~ 1.3. This is similar 
to the s cale height of the old ste llar population (z s =0.56 kpc in 
3.6 /mi, Holwe rda et alJlirTprepJ) in NGC 891. This is consistent 
with the results for low mass galaxies, where the vertical dust 
distribution has a sca le height somewh at larger than that of the 
main sequence stars dSeth et all 12005}) . However, it does com- 
plicate the picture for high mass galaxies, where previously it 
was thoug ht that the ISM comple tely collapsed into the typical 
dustlanes dDalcanton et alll2004 . 

The scale heights were calculated from profiles with a ver- 
tical resolution of l"to ensure that we had enough points to fit. 
These scale heights can be considered a lower limit to the real 
scale height of the dust since the scale height calculated in this 
way is a combination of the dust scale height and the Ha distri- 
bution. 

Higher up (|z| > 20") the profile hints at a second com- 
ponent on the East side (e.g., negative vertical offset) with an 
even larger scale height . This would imply that there is a 
dust halo with a scale height > 0.6 kpc. Thus, the extraplanar 
dust has a much larger scale height than previously assumed 
dXilouris et all [19991 iDalcanton et all l2004t ISethet all l200l 



Bian chUlin prep.1) . The quality of the data, in our opinion, is not 
good enough to fit a scale height to this second component. 

Moreover, if we fit the 24 /mi emission with a double ex- 
ponential we find the "normal" small scale height for the region 
close to the plane (z24/jm=0.22 ± 0.05 kpc) but a much larger sec- 
ond component (z24/im=E3 + 0.1 kpc) above 30"(1.4 kpc). These 
values are no more than indicative since we did not correct for 
the wings of the PSF in 24 //m image. 

The fact that the inferred values for r a dd are very reasonable, 
implies that the ionized gas in the halo is closely related to the 
star formation in the plane of the galaxy. Moreover, if our model 
is the correct view of NGC 891, the ionized gas and the attenuat- 
ing dust follow very different distributions in the halo: the dust is 
equally distributed around the disk while the ionized gas remains 
above the position it was initially ejected for a long time. 

5. Conclusions 

For the first time we quantitatively analyze the asymmetry 
between the North East and the South West side of the NGC 891 
at several wavelengths, thus tracing different components from 
the galaxy. We determine the ratio of asymmetry by dividing 
large strips on the North East side by the same area on the South 
West side of the galaxy (Fig. [2] upper panel). From this analy- 
sis it is immediately obvious that star formation tracers affected 
by dust attenuation (e.g., Ha, UV) show a much larger asymme- 
try than star formation tracers unaffected by dust attenuation (24 
/urn, radio continuum). We also see that the other components in 
the galaxy, like the old stellar population (3.6 /mi) and the neu- 
tral gas (HI) have only a small asymmetry. 

Assuming that 24 /mi emission from star forming regions 
dominates in the plane, t he 24 //m emission h as a direct correla- 
tion with star formation dCalzetti et all 120051) . The small asym- 
metry we find in the 24 /mi e mission is confirme d by the result 
found in the radio continuum dDahlem et all [19941) . which leads 
us to believe that 24 fim is indeed a good tracer of star forma- 
tion in NGC 891. Therefore, the asymmetry in Ha is most likely 
caused by dust in and above the plane. 

From this information we construct a simple symmetric 
model for NGC 891. This model is mainly based on the trail- 
ing spiral arm idea where the HII regions are located in front 
of the dust lanes (Fig [3]). We assume that the emission from the 
ionized gas is located in and above these HII regions but that the 
dust distribution in the halo is much more isotropic over the disk. 
This would lead to a longer sight line through the dust towards 
the Ha emission on the South West side of the galaxy. 

From this model we can derive the additional attenuation 
(i"add) in the halo between the spiral arms (Fig [2 lower panel). 
We see that at heights ~ 40" (1.9 kpc) the additional attenuation 
of Ha becomes negligible. This is either caused by the absence 
of dust at this height or by the distribution of the ionized gas 
becoming such that our lines-of-sight toward the ionized gas be- 
come equal on both sides. This last scenario, implies that the 
attenuation in the halo is even higher above the plane than de- 
rived from our simple model. 

The opacity values calculated in this letter only represent the 
additional opacity between the spiral arms (see fig. [3] straight ar- 
row. Thus, these values can be considered as a lower bound on 
the true opacity through the whole galaxy as a function of verti- 
cal height. 

Common wisdom still has it that extinction in the disk is pre- 
dominantly in the plane itself. Ty pical scale heights are much 
smaller than the stellar ones (Xilo uris etafl Il999t IBianchil 
lin prepj) or similar dSeth et ail 120051) . This paper shows that 



P. Kamphuis et al.: A dust component ~ 2 kpc above the plane in NGC 891 



5 



there is a second spatial component to the dusty ISM with a 
much more extended scale height. 

If our view of NGC 891 is correct then the unattenuated Ha 
emission above the plane is closely correlated with star form- 
ing regions in the plane. This implies that the ionized halo gas 
is mostly brought up from the plane through processes related 
to star formation. Furthermore the distribution of the attenuating 
dust above the plane is much more isotropic, and therefore more 
stable, than that of the ionized gas. This is probably caused by 
the ionized gas cooling down and recombining before it reaches 
this stable configuration. 
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